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A  Model  for  Predicting  Erosion  and 
Sediment  Yield  from  Secondary  Forest 
Road  Construction 

Charles  F.  Leaf1 

One  of  the  more  visible  and  controversial  environmental  impacts  associat- 
ed with  timber  harvesting  and  development  in  central  Colorado  is  road  con- 
struction. Better  tools  are  needed  to  quantify  the  effect  of  soil  disturbance  on 
erosion  onsite,  and  the  subsequent  yield  of  sediment  downstream. 

This  Note  summarizes  available  data,  and  from  this  base,  proposes  a  pre- 
liminary model  for  predicting  an  index  of  onsite  erosion  and  downstream 
sediment  yield. 
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Effect  of  Timber  Harvesting 
and  Road  Construction 

Watershed  erosion  and  damage  to  water 
quality  from  road  construction  and  timber  har- 
vesting can  be  significantly  reduced  through 
proper  planning,  construction,  and  followup 
maintenance  (Packer  and  Laycock  1969,  Meg- 
ahan  1972a).  For  example,  Leaf  (1970-71) 
showed  that  on  Fool  Creek,  in  the  Fraser  Experi- 
mental Forest,  road  construction  associated 
with  timber  harvesting  resulted  in  minimum 
erosion  damage  with  apparently  no  reduction  in 
water  quality.  The  3.3  miles  of  main  access  road 
were  carefully  located  to  avoid  the  stream  chan- 
nel and  to  minimize  erosion.  Timber  was  made 
accessible  by  an  additional  8.8  miles  of  spur 
roads  laid  out  along  contours.  Spur  roads  were 
provided  with  surface  drainage  and  culverts  at 
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Sediment  yield,  including  mineral  soil  and 
debris,  was  relatively  high  immediately  follow- 
ing road  construction.  However,  yield  de- 
creased rapidly  and  approached  preharvest 
levels  after  6  years  despite  a  persistent  increase 
in  runoff  caused  by  logging.  Other  studies  sup- 
port the  time  trend  results  observed  on  Fool 
Creek,  even  though  soil  and  geologic  conditions 
were  different  (Megahan  1974).  It  is  important 
to  understand  that  these  conclusions  apply  only 
to  surface  erosion,  and  not  to  mass  erosion, 
which  can  occur  through  the  disturbance  of 
very  steep  and  naturally  unstable  topography. 


Sediment  Yield  Model 

Megahan  (1974)  has  proposed  a  negative  ex- 
ponential equation  with  a  linear  component  con- 
taining three  parameters  to  describe  the  time 
trends  in  erosion  discussed  above.  The  equation 
can  be  expressed  as: 


E  =  ent  +  S0(l  -  e"kt) 


[1] 


where 

E  =  the  cumulative  onsite  erosion  at  time  (t) 
after  disturbance, 

en  =  an  estimate  of  the  long-term  "normal" 
erosion  rate  on  the  disturbed  area, 

S0  =  an  index  of  the  total  amount  of  soil  avail- 
able for  erosion  due  to  disturbance,  and 

k  =  an  index  of  the  rate  of  decline  of  erosion 
following  disturbance. 


According  to  Megahan,  (en)  is  the  long-term 
erosion  "norm"  which  is  reestablished  after  a 
site  is  disturbed.  In  some  areas,  this  new  norm 
may  be  higher  than  the  long-term  erosion  under 
natural  conditions  because  of  irreversible 
changes  in  site  factors.  However,  the  limited 
data  from  Fool  Creek  and  from  undisturbed 
watersheds  in  the  area  indicate  that  the  re- 
established erosion  norm  on  Fool  Creek  is  es- 
sentially the  same  as  before  disturbance.  There- 
fore, equation  [1]  is  used  to  predict  erosion 
time  trends — with  one  modification.  In  this  re- 
port, (en)  is  defined  as  the  natural  long-term  sedi- 
ment yield. 

The  logging  operation  on  Fool  Creek  required 
disturbance  of  35  acres  (Goodell  1958);  virtually 
all  of  this  involved  construction  of  roads  and 
landings.  Because  most  of  the  disturbed  area  is 
occupied  by  roads,  it  can  easily  be  described  in 
terms  of  road-design  variables  (fig.  1).  For  ex- 
ample, 12.1  miles  of  road  were  constructed 
across  slopes  which  averaged  26  percent.  Thus, 


the  35  acres  of  road  area  corresponds  to  an  "ef- 
fective" cross  section  having  the  following  char- 
acteristics: 


"effective  width" 
cut  and  fill  slopes 
width  of  disturbed  area  = 
area  disturbed  per  mile : 


14  ft 

1V2:1  (33.7°) 
23  ft 
2.8  acres 


These  results  were  obtained  from  the  following 
equations  which  can  be  derived  from  figure  1: 


D  =  W+-W/2  tanp 
tan 9F -  tanp 

when  9F  =  8C  =  0, 

tan  p 


D  =  W 


1  + 


W/2  tan  p 
tan6c  -  tanp 

D  ^  W;  6  2=  p 


[2] 


[3] 


tan  0  -  tan  p 

where 

D  =  total  disturbed  length  perpendicular  to  the 

centerline  in  ft, 
p  =  steepness  of  the  sideslope  in  degrees,  and 
0  =  angle  of  cut  and  fill  slopes  in  degrees 

(0C  =  angle  of  cut  and  0F  =  angle  of  fill). 

Equation  [3]  assumes  balanced  cut  and  fill 
(that  is,  that  the  centerline  bisects  the  road 
width).  This  is  not  usually  the  case,  since  the 
cross  section  can  vary  from  total  cut  to  total  fill 
in  actual  practice.  It  is  assumed,  however,  that  a 
sufficiently  accurate  index  of  the  total  area  dis- 
turbed can  be  obtained  by  estimating  an  "effec- 
tive" width  and  average  cut  and  fill  slopes  for 
the  proposed  road  system.  Such  estimates  re- 
quire considerable  judgment  and  a  knowledge 
of  the  topography. 

Equation  [3]  provides  a  means  for  expressing 
the  area  disturbed  in  terms  of  watershed  and 
engineering  design  parameters.  When  this  ex- 
pression is  used  in  combination  with  equation 
[1]  (provided  that  data  are  available),  it  should 


Figure  1 . — Simplified  cross  section,  illustrating  design 
variables. 
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give  the  land  use  planner  some  latitude,  subject 
to  the  limitations  and  assumptions  discussed 
below,  in  predicting  the  probable  erosion 
hazard  for  alternative  road  designs. 


Onsite  Erosion  Versus 
Sediment  Yield  Downstream 

To  fit  parameters  to  equation  [1],  the  sediment 
yields  measured  at  the  gage  on  Fool  Creek  were 
adjusted  to  unit  yields  per  acre  of  the  area 
actually  disturbed  (35  acres).  These  data  (sum- 
marized in  table  1)  represent  gross  volumes  of 
mineral  soil  and  debris  immediately  following 
road  construction.  It  should  be  emphasized  that 
table  1  summarizes  erosion  indices  and  not  the 
actual  erosion  which  took  place  on-site,  since 
the  data  were  collected  in  a  debris  pond  at  a 
point  downstream.  The  accumulated  erosion 
data  were  plotted  in  figure  2,  and  [1]  was  fitted 
using  a  nonlinear  procedure: 


Table  1 . --Summary  of  sediment  yields  from  Fool 
Creek,  Fraser  Experimental  Forest 


Year 

Sed i men  t 
yield 

Year 

Sed i ment 

•   i  j  l 
yield 

Year 

Sed  i ment 

■   i  j  l 
yield 

Ft3 /acre 

Ft  /acre 

Ft  /acre 

1952 

44.9 

1958 

42.8 

1964 

(3) 

1953 

22.6 

1959 

8.6 

1965 

25-5 

1954 

(2  3) 

I960 

13.9 

1966 

5.6 

1955 

(2  3) 

1961 

6.1 

1967 

6.4 

1956 

36.7 

1962 

16.3 

1968 

11.1 

1957 

70.1 

1963 

(3) 

1969 

13.9 

1970 

15.2 

Based  on  35  acres  disturbed;  average  dry  unit 
weight  =  85  lb/ft3  (range : 1 6- 1 20  lb/ft3). 

2Timber  harvest. 

3  Neg 1 i  g  i  b 1 e  accumulation. 


E 
E 


0.28t  +  401.3(1  -  e" 


[4] 


=  an  index  of  the  cumulative  onsite  erosion 
in  ft3/acre  at  time  (t)  after  road  construc- 
tion on  Fool  Creek, 

=  0.28  ftVacre  which  was  determined  from 
(1)  average  erosion  on  Fool  Creek  after  the 
6th  year,  and  (2)  long-term  yields  from  two 
undisturbed  watersheds,  and 

=  the  number  of  years  after  the  initial  dis- 
turbance. 


Three  assumptions  were  made  in  order  to 
develop  the  model  for  predicting  erosion  and 
sediment  yields.  First,  it  was  assumed  that 
equation  [4]  provides  a  better  index  of  erosion 
than  equations  based  on  rainfall-derived  credi- 
bility indices.  Such  indices  do  not  predict  time 
trends,  and  furthermore,  do  not  account  for  the 
effects  of  snowmelt,  which  is  responsible  for 
much  of  the  sediment  yield  from  the  densely 
forested  subalpine  zone  in  Central  Colorado. 
The  second  assumption  was  that  onsite  erosion 
is  proportional  to  the  area  disturbed.  Finally, 
it  was  assumed  that  the  delivery  ratio  is  constant 
for  a  given  watershed  size,  regardless  of  the 
amount  of  area  disturbed.  These  assumptions 
involve  complex  interactions  between  the 
hydrology,  geology,  and  soils,  which  need  to  be 
verified  by  additional  study. 

Based  on  the  simplifying  assumptions  dis- 
cussed above,  equations  [2]  and  [4]  were 
combined  in  order  to  predict  an  index  of  the 
cumulative  onsite  erosion  (S),  as  a  function  of 
the  width  of  roadbed  (W),  average  watershed 
sideslope  (p),  and  angle  of  cut  and  fill  (0C  and 
6F). 
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Figure  2. — Accumulated  erosion  as  a  function  of  time 
on  Fool  Creek,  Fraser  Experimental  Forest. 
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Thus 

S  =  0.121  DEn 

Substitution  of  equation  [2]  for  (D)  yields 
tan  p  r  1 


S  =  0.121  WEn 
1 


1  + 

] 


tanGc  -  tanp 


[5] 


tan  6F-  tanpJ 
where 

S  =  an  index  of  the  cumulative  erosion  in  Ft3, 
E  =  given  by  equation  [4], 
t  =  the  number  of  years  after  initial  disturb- 
ance, and 

n  =  the  number  of  miles  of  road  construction. 

Because  equation  [5]  is  expressed  in  terms  of 
engineering  design  variables,  its  use  should 
provide  an  indication  of  the  probable  erosion 
impacts  of  alternative  road  systems.  The  yields 
expressed  on  a  watershed  basis  are  given  by  the 
equation 


S 
A 


[6] 


where 

A  =  the  area  of  the  watershed. 

Equation  [6]  is  valid,  provided  that  the  up- 
stream drainage  area  does  not  exceed  1  mi2.2 
Sediment  yields  at  downstream  points  would  be 
less,  since  delivery  ratios  are  inversely  related 
to  watershed  area. 


Discussion  and  Conclusions 

One  of  the  more  significant  results  from  sedi- 
ment yield  studies  in  mountain  watersheds  is 
that  most  of  the  erosion  impact  occurs  within  a 
few  years  after  disturbance.  This  time  factor 
should  be  considered  in  land  use  planning  from 
two  standpoints:  protection,  and  long-term 
effects  on  hydrologic  parameters  such  as  water 
quality. 

Equations  which  require  erodibility  indices 
based  on  rainfall  intensity  may  be  grossly  in 
error  when  applied  to  much  of  the  subalpine 
zone  in  central  Colorado,  where  much  of  the 
sediment  yield  results  from  melting  snow. 
Equation  [1],  proposed  by  Megahan  (1974), 
appears  more  appropriate.  By  describing  the 
disturbed  area  in  terms  of  watershed  slope  and 
engineering  design  parameters,  as  given  in 

Equation  [4]  is  based  on  data  collected  from  a  714-acre 
experimental  watershed,  and  therefore  applies  only  to  drainage 
areas  of  approximately  1  mi2  in  size. 


equation  [5],  the  land  use  planner  has  some 
flexibility,  within  the  limitations  discussed 
above,  in  evaluating  the  potential  impacts  of 
various  road  systems. 

Because  the  coefficients  in  equation  [4]  were 
developed  from  a  very  limited  amount  of  data, 
they  may  not  be  generally  applicable  through- 
out the  Rocky  Mountain  Region.3  Nevertheless, 
they  are  based  on  the  best  information  available, 
and  should  be  considered  as  tentative  estimates 
until  more  data  become  available. 

The  model  described  in  this  report  is  based  on 
data  obtained  from  a  carefully  constructed  road 
system  and  a  high  standard  of  followup  main- 
tenance. Any  application  of  the  model  should 
presume  similar  standards  of  construction  and 
maintenance. 
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3/f  Is  worth  noting  here  that  soils  on  Fooi  Creek  are  uniform 
and  derived  from  gneiss  and  schist  rocks  (Retzer  1962).  They 
are  deep,  coarse  textured,  and  capable  of  absorbing  virtually 
all  of  the  water  during  peak  snowmelt. 
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